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ULTRASONIC  TRANSIT  TIME  MEASUREMENTS 
USING  FAST  ANALOG  TO  DIGITAL  CONVERSION 
AND  A  FAST  FOURIER  TRANSFORM  BASED  ALGORITHM 


by 


S.J.  Rumble  and  J.G.  Sparrow 


SUMMARY 

This  paper  describes  the  application  of  a  fast  (lOOMHz)  analog 
to  digital  voltage  converter  and  a  discrete  fast  Fourier  transform 
based  algorithm  to  the  determination  of  ultrasonic  transit  times.  The 
precision  obtained  was  typically  better  than  1  nanosecond. 

Treuisit  times  over  a  range  of  frequencies  {2-7MHz}  have  been 
measured  using  a  broadband  transducer  excited  by  a  chopped  CW  RF  pulse 
The  frequencies  chosen  corresponded  to  those  available  from  the 
discrete  fast  Fourier  transform. 

Documented  listings  of  the  programs  used  are  provided. 


COMMONWEALTH  OF  AUSTRALIA  1985 


1. 


INTRODUCTION 


The  potential  of  ultrasonic  methods  to  determine  residual 
stresses  has  been  recognized  by  numerous  workers  (1,2^3).  In  recent  years 
the  availability  of  improved  instrumentation  has  led  to  renewed  effort 
(4,5)  to  develop  practical  techniques  of  ultrasonic  determination  of 
residual  stress.  The  inproved  instrumentation  includes,  commercially 
available  broadband  shear  wave  trcmsducers,  EMAT'S  (electromagnetic 
acoustic  transducers) ,  fast  analog  to  digital  voltage  converters  and 
associated  microprocessor  based  hard%rare  and  software.  The  acquisition  of 
digital  records  of  the  voltage  response  of  the  broadband  transducer-sanple 
system  and  subsequent  conputer  processing  enables  the  frequency  dependence 
of  ultrasonic  transit  times  to  be  determined.  The  digital  records  will  also 
help  in  the  development  of  the  theory  of  propagation  of  ultrasonic  waves  in 
stressed,  textured  materials,  as  they  will  enable  a  more  detailed  conparison 
of  the  theoretical  and  experimental  results  to  be  made.  In  this  paper, 
details  are  given  of  the  application  of  specific  instrumentation  and 
software  to  the  determination  of  ultrasonic  transit  times  at  various 
frequencies. 


2 .  EXPERIMENTAL 


A  block  diagram  of  the  experimental  system  used  is  shown 
in  Figure  1.  The  transducer  was  a  Panametrics  type  V157  SMHz  broadband  shear 
wave  transducer  with  an  active  PZT  element  of  0.125  inch  diameter.  The 
transducer  was  excited  with  a  chopped  CW  RF  pulse  generated  by  a  HP3314A 
function  generator.  The  treinsducer  voltage  response  was  amplified  and  then 
digitized  at  10  nanosecond  intervals  by  a  Data  Precision  D6000  digital 
waveform  analyzer  with  a  Model  620  lOOHHz  8bit  11.28V  range  plug-in.  The 
Model  620  plug-in  could  acquire  up  to  30000  points  with  a  10  nanosecond 
interval  between  points,  which  corresponds  to  an  elapsed  time  of  300 
microseconds.  Both  the  D6000  and  HP3314A  had  remote  computer  control 
capedaility  and  were  interfaced  to  a  HP86  desktop  computer  via  an  IEEE-488 
bus. 


The  D6000  has  extensive  programmable  record  handling  and 
data  manipulation  capabilities.  In  particular,  it  was  able  to  perform 
discrete  fast  Fourier  transforms  (DFFT)  on  any  part  of  a  digital  record. 

The  HP86  was  used  as  the  controlling  computer  on  the  IEEE-488  interface 
bus.  The  D6000  was  initialized  so  as  to  be  triggered  by  the  transducer 
I  exciting  pulse,  and  to  take  sufficient  samples  to  encompass  all  return 

echoes  with  an  adequate  signal-to-noise  ratio.  Typically,  for  a  10mm 
thick  2024T4  aluminium  alloy  sample,  the  shear  wave  transit  time  was  6400 
nanoseconds  and,  with  the  broadband  transducer  used  the  signal  to  noise  ratio 
was  falling  rapidly  past  echo  4.  An  example  of  a  digital  voltage  record 
for  a  10mm  path  in  2024T4  aluminium  alloy  and  at  4.69MHz  is  shown  in 
Figure  2. 


The  ultrasonic  transit  time  was  obtained  from  the  digital 
record  by  determining  the  corresponding  points  of  the  same  phase  (to 
within  an  integral  multiple  of  21)  in  the  return  echoes  of  interest.  In 
the  measurement  of  transit  time,  a  precision  of  0.1  to  1  nanoseconds  was 
required  to  determine  the  stress  level  to  IMPa.  The  matching  of  corresponding 
points  in  different  echoes  is  con^licated  by  the  frequency  dispersion  of 
the  ultrasonic  velocity  and  Papadakis  (6)  has  noted  that  care  needs  to  be 
taken  to  ensure  that  the  transit  time  corresponding  to  the  phase  velocity 
and  not  the  group  velocity  is  measured. 


In  the  system  described  here,  the  transit  time  was  determined 
by  initially  selecting  from  the  digital  records  the  small  regions  within 
the  record  which  contained  the  echoes.  The  number  of  sample  points,  and 
hence  the  time,  of  these  regions,  and  the  time  at  which  these  regions 
started  were  parameters  in  the  main  controlling  program  on  the  HP86.  These 
parameters  were  subsequently  sent  by  the  HP86  to  the  D600  as  part  of  the 
D6000*s  controlling  program. 


Using  the  terminology  introduced  by  Allen  et  al  (5) ,  the 
time  between  these  start  points  is  referred  to  as  the  "coarse"  transit 
time.  A  fine  transit  time  adjustment  was  calculated  by  first  performing  a 
DFFT  on  the  individual  echo  regions.  The  phase  for  the  appropriate 
frequency  calculated  by  the  DFFT  for  each  of  the  two  echoes  was  subtracted, 
and,  this  difference,  on  multiplication  by  the  angular  frequency  was  the 
fine  transit  time.  The  transit  time  was  then  the  sum  of  the  coarse  transit 
time  and  the  fine  transit  time. 


The  DFFT  required  that  the  number  of  points  be  a  power  of 
two.  Because  of  windowing  effects  the  region  was  chosen  so  as  to  encompass 
the  entire  echo.  This  constraint,  along  with  execution  time  constraints, 
put  a  limit  on  the  number  of  RF  cycles  in  the  chopped  CW  RF  pulse.  For  the 
frequency  range  possible  (2-7MHz)  and  for  a  256  point  DFFT,  5  or  6  RF  cycles 
could  be  used.  As  the  time  length  of  the  DFFT  record  determined  the 
frequency  resolution,  the  choice  of  a  256  point  DFFT  with  10  nanosecond 
sampling  time,  constrained  the  frequencies  that  could  be  used  to  integral 
multiples  of  0.391MHz. 


3.  DESCRIPTION  OF  COMPUTER  PROGRAM 


The  main  controlling  program  on  the  HP86  prompted  the 
user  to  enter  the  following  parameters:  the  length  of  the  DFFT  in  sample 
points,  the  start  points  of  all  the  regions  within  the  record,  and  the 
number  of  records  per  frequency  to  be  averaged.  The  program  determined 
the  appropriate  integer  for  the  start  and  finish  frequency,  given  the  DFFT 
length,  and  assuming  10  nanoseconds  sampling  time,  from  which  the  frequency 


resolution  could  be  calculated.  It  then  initialized  the  HP3314A  and 
initialized  and  programmed  the  D6000.  Vfhen  this  was  con^leted,  the 
D6000  was  requested  to  acquire  a  record  and  to  run  its  program,  which 
selects  and  performs  DFFT's  on  the  required  regions.  The  HP86  waited 
for  this  program  to  finish,  and  then  accessed  the  phase  information  for  each 
echo  at  the  appropriate  frequency,  and  checked  to  see  if  it  was  valid. 

The  phase  differences  were  calculated  and  converted  to  fine  transit  times, 
which  was  followed  by  the  updating  of  the  running  sums  for  the  calculation 
of  the  average  and  standard  deviation.  The  cycle  of  acquisition  of  data 
record,  DFFT's  by  the  D6000,  transfer  of  phase  data  to  HP86,  and 
calculation  of  fine  transit  times  for  each  echo  pair  was  repeated  the 
requested  number  of  times  for  each  frequency.  The  program  concluded  by 
printing  all  the  relevant  information.  A  documented  listing  of  the  HP86 
program  and  the  06000  program  is  contained  in  Appendix  1. 


4 .  RESULTS 


The  technique  detailed  above,  has  been  used  to  investigate 
the  frec[uency,  temperature  and  stress  dependence  of  the  ultrasonic  transit 
times  in  a  extruded  10mm  2024T4  aluminium  alloy  sample.  The  results  of 
these  investigations  are  reported  elsewhere  (7) .  An  example  of  the  output 
of  the  program  is  given  in  Appendix  1.  Examples  of  plots  of  the  actual 
data  in  a  region  encompassing  an  echo  and  overplots  of  sinusodial  cvurves  of 
the  phase  calculated  by  the  DFFT  and  with  the  frequency  of  the  exciting  pulse 
are  shown  in  Figures  3,4,  and  5. 


5.  DISCUSSION 


The  standard  deviation  of  the  mean  of  the  transit  time 
measurements  have  typically  been  below  1  nanosecond  for  echo  1  to  2  and 
echo  2  to  3.  This  precision  allows  determination  of  stresses  to  the  order 
of  IMPa.  However,  for  later  echo  pairs,  the  standard  deviations  can  increase 
to  above  1  nanosecond.  At  this  stage  the  transit  times  are  not  corrected 
for  diffraction  effects  (6,8),  and  any  211  errors  (8)  are  corrected  by  hand. 
The  effect  of  diffraction  is  to  produce  different  transit  times  between 
echoes  1  and  2  to  those  between  echoes  2  and  3  etc..  The  diffraction 
corrections  calculated  by  Allen  show  that  the  transit  times  between  the 
early  echoes  are  too  short,  but  the  correction  quickly  reduces  to  zero  for 
later  echo  pairs.  However,  even  after  the  addition  of  these  corrections, 
Allen  (8)  has  observed  that  the  transit  times  between  the  various  echo  pairs 
are  often  different. 


In  the  results  obtained  using  the  above  technique,  there  was 
a  significant  variation  in  the  transit  times  between  the  various  echo  pairs, 
which  could  not  be  explained  by  the  diffraction  effect  calculated  by  Allen (8) 


There  are  a  number  of  other  factors  which  could  cause  these  transit  time 
variations.  These  include,  interference  from  shear  to  longitudinal 
mode-converted-waves,  interference  due  to  echoes  from  side  walls,  generation 
of  incompletely  linearly  polarized  shear  waves  by  the  transducer,  complex 
interactions  arising  from  the  texture  in  the  sample  emd  the  transducer/ 
sample  bond,  and  incomplete  dancing  within  the  backing  of  the  broadband 
transducer.  The  multiplicity  and  conqplexity  of  the  possible  factors 
involved,  highlights  the  advantages  of  having  a  complete  digital  record 
of  the  transducer /sasple  system.  This  is  because  it  will  enable  computer 
optimization  of  the  various  parameters  of  amy  theoretical  model  by 
comparison  with  the  experimental  detezmined  digital  records.  Further 
experimental  and  theoretical  study  will  be  carried  out  to  attempt  to 
resolve  the  above  problems.  Currently  experimental  work  is  being  performed 
with  the  aim  of  completely  utilizing  the  broadband  nature  of  the  transducer. 
In  this  system,  the  transducer  is  excited  by  a  very  fast  rise  time  pulse. 

The  generated  ultrasonic  pulses  have  significant  energy  over  a  broad  range 
(2-7MHz)  of  frequencies,  and  after  Fourier  analysis,  the  phase  of  the 
various  frequency  conqponents  can  be  determined.  In  this  way,  the  transit 
times  at  a  number  of  frequencies  can  be  determined  simultaneously.  The 
results  of  this  %«ork  will  be  reported  in  due  course. 


6.  CONCLUSION 


This  work  has  demonstrated  the  potential  of  digital 
analysis  techniques  in  the  ultrasonic  determination  of  residual  stresses. 

The  method  described  here  has  enabled  ultrasonic  transit  times  to  be 
determined  to  a  precision  of  typically  1  nanosecond  over  a  range  of 
frequencies.  The  transit  time  measurements  obtained  with  this  method  have 
shown  that  further  theoretical  and  experimental  work  is  required  before  the 
effect  of  diffraction,  mode  conversion,  texture  and  frequency,  among  others, 
are  fully  understood. 
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REGION  FROM  FIGURE  2  USED  FOR  ECHO  3  AND  AN  OVERPLOT  OF  SINUSOIDAL  CURVE 
THE  PHASE  CALCULATED  BY  THE  DFFT. 


APPENDIX  1 


Listing  of  the  main  controlling  program  on  the  HP86 
the  program  on  the  06000  and  typical  output  from  the  RP86  program. 
The  output  is  for  a  10mm  transit  path  in  2024T4  aluminium  alloy.  N< 
that  the  transit  times  are  not  corrected  for  2n  errors,  nie  date 
in  the  fourth  column  in  the  output  are  calculations  of  the  standard 
error  of  the  mean. 

HP86  program  name  :  TTIME. 


10  DIM  PH(5.20) ,FINETM<5,20) .ES(5) .DIFFCS) .TIM(5) 

20  REM  PH(I.J)  CONTAINS  THE  PHASE  FOR  ECHO  I.  RECORD  J 

30  REM  FINETTd.J)  CONTAINS  THE  FINE  TRANSIT  TIME  BETWEEN  ECHO  I  +  l  AND 

40  REM  ECHO  I  FOR  RECORD  J 

50  REM  ESd)  CONTAINS  START  IN  DATA  POINTS  OF  Ith  ECHO  REGIONS 

60  REM  IN  DATA  RECORD  ON  D6000 

70  REM  DIFFd)  CONTAINS  PHASE  DIFFERENCE  FOR  ECHO  PAIR  I  +  l  AND  I 

80  REM  TIMd)  CONTAINS  TOTAL  TRANSIT  TIME  FOR  ECHO  PAIR  I-H  AND  1 

90  DIM  DAT»C53,ANGLE*C53.TEMP*C53.STRESS*t63 

100  REM  STRING  VARIABLES  USED  TO  STORE  DATE,  ANGLE  OF  TRANSDUCER, TEMPERATURE 
no  REM  AND  STRESS 

120  STARTF=2.73  '  STARTING  FREQUENCY  IN  MHz 

130  ENDF=5.86  !  FINAL  FREQUENCY  IN  MHz 

140  DISP  "ENTER  LENGTH  OF  FFT  (POWER  OF  TWO)  " 

150  INPUT  LENGTH 

160  DELF= 100 /LENGTH  •  FREQ  RESOLUTION  IN  MHz  ASSUMING  lO  NANOSECS 
170  NST=FNRD (STAPTF/DELF,0)  !  LOCATION  OF  DATA  FOR  STARTF  IN  FFT  OUTPUT 
ISO  NEND=FNRD (ENDF/DELF.O)  '  LOCATION  OF  DATA  FOR  ENDF  IN  FFT  OUTPUT 
190  DISP  "INPUT  NUMBER  OF  RECORDS  PER  FREQ  FOR  AVERAGING" 

200  INPUT  JMAX 

210  ERPNUM=JMAX/2  '  NUMBER  OF  INVALID  RECORDS  ACCEPTED  BEFORE  USER  IS 

220  REM  PROMPTED 

230  PEM  PROMPT  FOR  RUN  INFORMATION 

240  DISP  "ENTER  DATE  .  ANGLE  AND  TEMP, STRESS" 

250  INPUT  DAT*, ANGLE*. TEMP*, STRESS* 

260  PRINT  DAT*S<"  CALCULATION  OF  MEAN  AND  S.D  OF  FINE  TRANSIT  TIMES  FOR";JMAX;"  F 
ECORDS" 

270  PRINT  "USING  FUNCTION  GENERATOR— ROTATION  OF  TRANSDUCER  "S<ANGLE*»-"  DEGREES" 
280  PRINT  "SAMPLE  TEMPERATURE  "«<TEMP*.»-"  DEGREES" 

290  PRINT  "STRESS  ".STRESS* 

300  PRINT  C*'  PRINT  "ECHO#"  .  "FREQ  (MHz  )".  "MEAN  (NANOSECS)  ",  "STD  DEV/ (  "  s  JMAX  s  "  )  '  .  5" 
310  1=1 

^20  REM  PROMPT  FOR  STARTING  POINTS  OF  ECHO  REGIONS  WITHIN  DATA  RECORD 
330  DISP  "INPUT  START  OF  ECHO#": I:"  (-1  TO  FINISH)" 

340  I NPUT  ESd' 

350  IF  ESd)<  :•  -1  THEN  I-I  +  l  R  GOTO  330 
360  NUMECH0=I-1  '  NUMBER  OF  ECHOES  ENTERED 

370  REM  CHECF  TO  SEE  IF  D6000  REQUIRES  INITIALIZATION  AND  PROGRAMMING 
380  DISP  "ENTER  1  IF  DP6000  ALREADY  SETUP" 

390  INPUT  SETFLAG 

400  IF  SETFLAG=1  THEN  610 

410  NL=2'"’ 

420  FDR  1=1  TO  NUMECHO 

430  REM  INITIALIZE  DATA  RECORDS  ON  D6000 

440  CALL.  "MALE"  <  "ECHO"«.VAL»  d  ). LENGTH .  O ..  0000000 1  .  1 . 0 .  1 . 28 , 0  ) 

450  REM  SEND  PROGRAM  TO  D6000 

460  PEM  SELECT  ECHO  REGIONS  FROM  DATA  RECORD 

470  CALL  "SEND"  <  VAL*  d*NL)?"  ECHO"*-VAL*  d  )  ?>"=BUF .  A1  (  "?-VAL*  (ES  d  "«-VAL* 

(ES  d  '  *255'  !•'"  )  "  ) 

480  REM  PERFORM  DFFT  ON  ECHO  REGION 

490  CALL  "SEND"  <  VAL*  d*NL  +  IO'?-"  MAGS  "?'VAL  *  d  )  S- "  C  RHAS  " -VVAL  *  d  )  !•"  I=FFT  (ECHO"? 


VAL  *  d  '  .  6 . 0 . 0 .  1  2 '  "  ' 

500  next  I 

510  REM  INITIALIZE  HF3314A  FUNCTION  GENERATOR 

520  OUTPUT  307  ; "RC5"  •  RESTORES  HP3314  TO  STATE  STORED  IN  MEMORY  5 

530  REM  'MODE=  N  CYCLE 

540  REM  'AMPLITUDE  =  10.00  VOLTS 

550  REM  'Nl.lMBER  OF  CYCLES  =  5 

560  PEM  'TRIGGER  =INTEPNAL 

570  pem  INITIAIIZE  D60i:iO 

5po  roll  "SEND"  '  "FFnMFT=3:EPPM=3; TFGM= 1 ; DAPM: E>FroN= I : "  ' 


Cftl.L  "SEND"  <  "DSPM=1;  )'0FF*6.0li5;TRGLEV=1  .l''''.';NFTS' 1  »=409to:FEF:I0D' 1  >  =  l'-'nS"  ' 

600  FEM  LDDF  FOR  EACH  FPEOUENCv  FROM  N5T  TO  NEND 
610  FOP  N=N5T  TO  NEND 
620  REH  LOOP  FOP  EACH  ECHO  REGION 
630  FOP  1-1  TO  NUMECHO 

640  PEM  INITIALIZE  VARIABLES  USED  FOR  MEAN. STD.  DEV.  AND  ERROR  COUNTS 
650  6-  X2<I>“0  $  ERRCNT<I)*0 

660  NEXT  I 

670  FREO=DELF*N  ‘  CALCULATE  ACTUAL  FREQUENCY  IN  MHr 
680  REM  ROUND  FREQ  TO  TWO  DECIMAL  PLACES 
690  FREO=INT  (FREQ*IOO*.5) /lOO 

700  REM  SET  FUNCTION  GENERATOR  TO  REQUIRED  FREQUENCY 
710  OUTPUT  307  ;  •■FRe"«-VAL*  <FREO>  J-"MZ" 

•X20  REM  LOOP  FOR  THE  REQUIRED  NUMBER  OF  RECORDS  TO  BE  USED  FOR  AVERAGING 
730  FOP  J»0  TO  JMAX-1 

740  REM  ACQUIRE  DATA  RECORD  AND  RUN  PROGRAM  ON  D6000 
750  CALL  "SEND"  (  "ARM"  ) 

760  CALL  "SEND"  <  "DISARM"  ) 

770  CALL  "SEND"  (  "RUN"  » 

780  REM  LOOP  TO  SEE  IF  PROGRAM  HAS  FINISHED 
790  CALL  "GETVAL"  (  "PGMST" . PGMST  ) 

800  IF  PGMST =2  THEN  790 

810  REM  LOOP  FOR  EACH  ECHO  REGION 

820  FOP  I-l  TO  NUMECHO 

830  PEM  GET  Nth  POINT  IN  RECORD  "PHAS".  WHICH  CONTAINS  THE  PHASE 
840  CALL  "SEND"  (  "M*PHAS"«-VAL»  <  I )  ?• "  '  "*.VAL  *  (N>?")"  ) 

850  CALL  "GETVAL"  (  "M",M  > 

860  IF  M'O  THEN  M=360-M  8  GOTO  890 

870  M=-M  '  CONVERT  TO  POSITIVE  DEGREES  TO  FIRST  PEAK  PAST  T*0 
880  REM  STORE  PHASE  DATA 
890  PH(I.J'=M 
900  NEXT  I 

910  PEM  LOOP  FOP  EACH  ECHO  PAIR 
920  FOP  1=1  TO  NUMECHO- 1 
930  REM  COMPUTE  DIFFERENCE  IN  PHASE 
940  DIFF ( I t  =PH (1+ 1 . J I -PH< I , J  > 

950  REM  CHECK  TO  SEE  IF  RECORD  WAS  VALID 

960  PEM  ONLY  CHECK  DIFFERENCE  BETWEEN  ECHO  1  AND  ECHO  2 
970  IF  I  •=  2  then  1 180 

980  PEM  STORE  DIFFERENCE  IF  IT  IS  THE  FIRST  RECORD 
990  IF  J-0  THEN  D IFFCK  < I » -DIFF i I > 

1000  PEM  ACCEPT  RECORD  IF  VARIATION  FROM  FIRST  RECORD  DIFFERENCE  IS  OP  =  40 
1010  IF  ABS  'DIFF (I » -DIFFCK  ( I )  '  =  40  THEN  1 1 B- • 

1020  PEM  ACCEPT  RECORD  IF  AFTER  A  2P1  ADJUSTMENT.  THE  VARIATION  FROM  THE 
1030  PEM  FIRST  RECORD  DIFFERENCE  IS  ■  OR  «  2'> 

1040  IF  APS  'DIFF ' I  I +36"-DIFFCK  ' I >  '  «  2"  THEN  RH ■  I ♦ )  . J ' =PH ' I ♦ 1 . J ' *360  8  DIFFiI'  = 
DIFF(I'*360  8  GOTO  IIB'" 

1M5.-1  IF  ABS  'DIFF  '  I  ' -36-:— DIFFCk  ■  I  )  1  «  D-:-  THEN  PH  -  I  ♦  1  .  J  '  »PH  i  I  ♦  I  .  J  ' -36'.'  8  DIFF'I’  = 
DIFF 'I' -360  8  GOTO  1  1  e-- 

106'-'  REM  RECORD  REJECTED  INCREASE  ERROR  COUN  .  Bv  ONE 
1 07':'  EPRr NT  -  I  '  *ERRCNT  •  I  >  ♦  1 

1''80  REM  IF  ERROR  COUNT  IS  LESS  THAN  LIMIT  THEN  OBTAIN  A  NEW  RECORD 
1090  IF  EPRCNT  '  I  1  ■  «  ERRNUM  then  7«..;. 

11, ->0  REM  ERROR  COUNT  GREATER  THAN  LIMIT.  INFORM  USER 

1110  DISR  "MORE  THAN  ";ERPNUM;"  RECORDS  OUT  OF  ";J:"  RECORDS  HAVE  HAD  PHASE  DIFF 
S  " 

112''’  DISR  "GREATER  THAN  10  DEGREES  Frqm  the  INITIA  PHASE  DIFF" 

1130  DISR  "BETWEEN  ECHOES  ":I;"  AND  ";I*I 

1140  DISR  "ENTER  2  TO  STOP  .1  TO  INFREASE  ERROR  CHECK R.O  TQ  GET  new  CHECK  DIFF" 
115':'  INPUT  TFLAG8  IF  TFLAG-C  THEN  STOP 
1  1  v:'  IF  TFLAG-'"'  then  63  " 

1170  DISR  "ENTER  NEW  ERROR  CHECK  NLIMBEP"  8  INPUT  ERRNUM 
1  18':-  NEXT  I 

11®';'  REM  UPDATE  P'lNNINC-  SUMS  QF  PHASE  DIFFS  AMD  PHASE  DIFPS  SQUARED  FOR 
lO'-"-'  REM  AfrcERTED  RECORDS 
1  -  1  pnc  T  =  i  Tn  Ni'NPrwn-, 


122>.'  »  1  I  J  =y  '  n  +D1FF  '  I  I 

1230  X2'I>=y2<I)+DIFF(I>»DIFF<I) 

1240  NEXT  I 
1250  REM 

1260  NEXT  J  •  LOOP-  FOR  NEXT  RECORD 
1270  REM 

1280  REM  CALCULATE  MEAN  AND  STANDARD  DEVIATION  OF  MEAN.  AND  CONVERT  DEGREES 

1290  REM  TO  NANOSECS 

1300  FOR  1=1  TO  NUMECHO-1 

1310  X (I ) =X ( I) /JMAX 

1320  X<I)=X(I)/<.36»FREQ) 

1330  X2( I )=X2( I ) / ( . 1296*FRE0*FREQ) 

1340  SD<I)=X2<I)-JMAX»xa)*X<I) 

1350  SD  ( I  )  =  (SD(  I )  /  <  JMAX-1 )  l-'.S 

1360  PRINT  8  PRINT 

1370  REM  ROUND  VALUES  AND  PRINT 

1380  X=FNRD(X  <  I  )  ,  1  )  8  S=FNRD  <SD  <  I )  /  JMAX''.  5. 2) 

1390  PRINT  I;"  AND  “ ; I  + 1 . FREQ,  X  . S 
1400  FINETM<I,N)=X<I) 

1410  NEXT  I 
1420  REM 

1430  NEXT  N  •  LOOP  FOR  NEXT  FREQUENCY 
1440  REM 

1450  REM  PRINT  FINAL  RESULTS 
1460  PRINT 
1470  PRINT 

1480  PRINT  "OVERALL  TRANSIT  TIMES  BETWEEN  ECHOES  AND  DIFFERENCES" 

1490  PRINT  8  PRINT  "ECHOES" , "FREQ/ /PERIOD <nS >". "TRANSI T  T IMES (nS) ". "DIFFERENCES 
1500  FOR  N=NST  TO  NEND 
1510  DELF= 100 /LENGTH 
1520  FOR  1=1  TO  NUMECHO-1 

1530  REM  COMPUTE  TRANSIT  TIMES  FOR  EACH  ECHO  PAIR 
1540  TIM(I>=<ES<I+1)-ES(I) > * 10+FINETM < I , N) 

1550  NEXT  I 

1560  FOP  1=1  TO  NUMECHO-1 
1570  FREQ=N*DELF 
1580  S*="  " 

1590  IF  I<NUMECH0-1  THEN  S*=VAL*  <FNRD ' T IM ( 1+ 1 ) -T IM ( I ) . 1 ) ) 

1 600  PEP=FNRD ' 1 000 /FREQ. 1 > 

1610  REM  POUND  AND  PRINT 

1620  FRE0=FNRD(FREQ.2)  8  T=FNRD < TIM  a ) .  1 ) 

1630  PRINT  I; "AND"; I+1,FREQ.T 
1640  PRINT  "  ".PER."  " . S* 

1650  NEXT  I 
1660  NEXT  N 

1670  REM  PROMPT  USER  FOP  NEW  PARAMETERS  IF  REOUIRED 
1680  CISP  "ANOTHER  RUN  (1.0)" 

1690  INPUT  FLG 

1700  IF  FLG=0  THEN  STOP 

1^10  DISP  "ENTER  NEW  ANGLE. NEW  TEMP. NEW  STRESS" 

1720  INPUT  ANGLE*. TEMP*. STRESS* 

1730  DISP  "AD3UST  PRINTER  PRESS  < CONTV  TO  CONTINUE" 

1740  PAUSE 
1750  GOTO  260 
1760  STOP 
177i:i  REM 
1780  REM 

1760  REM  SUBROUTINE  TO  POUND  PARAMETER  IN  "VALUE"  TO  "PLACES"  PLACES 
1800  DEF  FNRD 'VALUE. PLACES' 

1810  VALUE=INT  (VALUE* 10  PLACES*. 5i /lO  PLACES 
1820  FNRD= value 
1830  FN  END 
1840  END 


10  SUB  "GETVAL"  (NAME*. VALUE) 

20  REM  GETVAL  OBTAINS  THE  VALUE  OF  THE  PARAMETER  ON  THE  D6000 
30  REM  WHOSE  NAME  IS  CONTAINED  IN  THE  STRING  NAME* 

40  OUTPUT  31 S  ;NAME* 

50  ENTER  315  ;  ANSWER* 

60  ENTER  315  ;  B* 

70  VALUE=VAL  (ANSWER*) 

BO  SUBEXIT 
90  SUBEND 


10  SUB  "HAKE"  <NAME*. LENGTH, XOFFSET.XPERPNT,XUNIT,YOFFSET.YFSD,YUNIT) 

20  REH  CREATES  A  DATA  RECORD  ON  THE  D6000 

30  REH  CONSULT  D6000  USERS  MANUAL  FDR  EXPLANATION  OF  PARAHETERS 

40  CALL  "SEND"  <  "PR0HPT=3; ERRM=3"  )  !  INITIALIZES  PPOHPT  AND  ERROR  HODE 

50  CALL  "SEND"  (  "TRCSRC ( 1 ) »BUF. A1 "  )  !  SETS  TRACE  ON  SCREEN  TO  BUF.Al 

60  OUTPUT  315  5"NPTS<1)"  '  PROHPTS  D6000  TO  SEND  NUHBER  OF  POINTS 

70  ENTER  315  ;  B»  '  IN  DATA  RECORD  BUF.Al 

80  LTH=OAL  (B*i  '  NUHBER  STORED  IN  LTH 

90  ENTER  315  ;  B*  '  RECEIVE  PROHPT 

100  CALL  "SEND"  <  "NPTS ' 1 > ="*«VAL*  (LENGTH)  )  !  ALTER  LENGTH  OF  BUF.Al 

no  XOFF=XOFFSET/ (LENGTH* XPERPNT)  •  CALCULATE  VALUE  TO  BE  SEND  TO  D6000  FOR 
120  REM  AN  OFFSET  OF  XOFFSET 

130  CALL  "SEND"  (  "TEMP1=BUF. A1 "  )  '  CREATE  A  TEMPORARY  FILE  THE  SAME  SIZE  AS 

140  REM  BUF.Al 

150  CALL  "SEND"  (  "TEMPI  (O.  “«-VAL*  (LENGTH-1  )*'")  =0. 0"  )  '  ZERO  ALL  VALUES 

160  CALL  "SEND"  (  "TEMP2=UNIT  (TEMPI .  O,  "fcVAL*  (  XOFF )  8<"  .  "?-VAL*  ( XPERPNT  >  S."  ,  "?-VAl  * 

(XUNIT)?.")  "  ) 

170  REM  ALTER  X  OFFSET  AND  X  UNITS  PER  POINT 

180  CALL  "SEND"  (  NAME**<"»UNIT  (TEMP2,  1 .  "*<VAL*  ( YOFFSET)  .  "?/VAL»  ( YFSD)  .  "J'VAt  * 
(YUNIT)*-")"  ) 

190  REM  RENAME  RECORD  AND  ALTER  Y  OFFSET  AND  Y  F.S.D 
200  CALL  "SEND"  (  "DEL  TEMPI ;DEL  TEMP2"  ) 

210  REM  REMOVE  TEMPORARY  FILES  CREATED 
220  CALL  "SEND"  (  "NPTS  ( 1 )  *"*<VAL*  (LTH)  ) 

230  REM  RESTORE  LENGTH  OF  BUF.Al 
240  CALL  "SEND"  (  "LOCAL"  ) 

250  REM  RESTORE  LOCAL  CONTROL  OF  D6000 
260  SUBEXIT 
270  SLIBEND 


10  SUB  "SEND"  (ft») 

20  REM  SENDS  STRING  IN  TO  D6000 

30  REM  SEND  ASSUMES  THAT  D6000  HAS  BEEN  INITIALIZED  TO  SEND  A  "  " 

40  REM  PROMPT  ON  ACCEPTANCE  OF  MESSAGE 

50  REM  INITIALIZATION  IS  PERFORMED  BY  CALL  "SEND"  ( "PROMPT-T" ) 

60  DISP  A* 

70  REM  DISPLAY  STRING  SENT  ON  HPB6  SCREEN 
80  OUTPUT  315  ; A» 

90  REM  315  =  DEFAULT  ADDRESS  OF  D6000  ON  IEEE-488  INTERFACE  BUS 
100  ENTER  315  ;  B» 

110  IF  B*m  =  ">"  THEN  SUBEXIT 
120  DISP  B« 

130  REM  DISPLAYS  ON  HP86  SCREEN  ANY  MESSAGE  SENT  BY  D6000  IN  RESPONSE  TO  STRING 
140  GOTO  100 
150  SUBEND 
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